A Single-Ion Nuclear Clock for Metrology at the 19 Decimal Place 
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The 7.6(5) eV nuclear magnetic-dipole transition in a sing le 229 Th 3+ ion may provide the founda- 
tion for an optical clock of superb accuracy. A virtual clock transition composed of stretched states 
within the 5F 5 / 2 electronic ground level of both nuclear ground and isomeric manifolds is proposed. 
It is shown to offer unprecedented systematic shift suppression, allowing for clock performance with 
a total fractional inaccuracy approaching lxlO -19 . 



The development of optical atomic clocks has produced 
tremendous advances in time and frequency metrology 
and their applications in tests of fundamental physics. 
Various ionic and neutral atomic species are now utilized 
for frequency standards with fractional inaccuracies at 
or below 10 -15 [TH5]. The most accurate systems ap- 
proach and even surpass the 10 -17 mark [TTJflT]. How- 
ever, the difficulty in simultaneous minimization of ex- 
ternal field and motional shift uncertainties, as well as 
clock linewidth limitations, make the next two orders of 
magnitude of inaccuracy reduction challenging. In 2003, 
Peik and Tamm pointed out that many of the standard 
external field shifts in atomic clocks could be suppressed 
or avoided in a 229 Th 3+ optical nuclear transition by uti- 
lizing the 2 Si/2 electronic level within both the nuclear 
ground and isomer manifolds |12j . This pioneering pro- 
posal illuminated a way towards a next-generation clock, 
however, their system is still subject to two important 
ion-clock limitations: the proposed mj? = 0f>mf = 
clock transition experiences a significant second-order 
differential Zeeman shift (~ 70kHz/mT 2 ), while electric- 
quadrupole decay of both clock states limits the clock 
linewidth to ~1 Hz. 

In this Letter, we show that the pair of stretched 
hyperfine states within the 5F 5 / 2 electronic ground 
level of both nuclear ground and isomeric manifolds in 
229 Th 3+ provide large suppression of all external field 
clock shifts, including the quadratic Zeeman shift, while 
providing a clock linewidth of <100/iHz. The clock 
is composed of the \5F 5 / 2 , I g = 5/2; F = 5, rap = ±5) 
O \5F 5/2 , I m = 3/2; F = 4, m F =±4) transitions (Fig. [T]) 
which, when averaged within a weak magnetic 
field, form a virtual clock transition. In gen- 
eral, stretched eigenstates in the coupled basis 
are also eigenstates in the uncoupled basis, i.e., 
\F = I+J ) m F =±F) = \J,±J)®\I,±I). For a struc- 
tureless point-like nucleus, both ground and excited clock 
state shifts are then exactly equal. When nuclear size 
and multipole moment differences in the clock levels are 
taken into account, differential shifts do arise, however 
they are suppressed well below levels of optical atomic 
clocks. This suppression encompasses shifts due to mag- 



netic fields, electric fields, and electric field gradients, 
as described in detail below. With such resilience to 
external-field perturbations, our proposal opens the door 
for metrology at a level unattainable by any other cur- 
rently proposed clock. 

Because the 5F 5 / 2 orbitals are the electronic ground 
levels, clock interrogation time is fundamentally limited 
only by the 10 3 - 10 4 s nuclear isomer excited state life- 
time [13 15 , providing potential for a resonator qual- 
ity factor of 10 19 - 10 20 . In a system where quantum 
projection noise dominates clock operation, the Allan 
deviation is a y {r) w ("fcik/^cik){T C ycie/T) 1/2 p2], where 
7 c ;fc is the clock excitation linewidth, is the clock 
transition frequency, T cyc i e is the clock cycle time, and 
r is the total integration time. With an ultra-stable 
interrogation laser and efficient state detection allow- 
ing for 7 cife fa 1/T cyc i e ps 1 Hz, a 229 Th 3 + clock at 1.8 PHz 
could produce a quantum-limited fractional instability 
of lxl0~ 15 t -1 / 2 , allowing for a statistical uncertainty 
of lxlO -18 to be achieved in 2 weeks of averaging. This 
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FIG. 1. Partial energy-level diagram of the 5F$/2 electronic 
ground levels within the nuclear ground and isomer manifold 
of 229 T n 3 +. The nuclear ground state hyperfine structure is 
taken from Ref.|17j. whereas the isomer hyperfine splittings 
are estimates. See text for details. 
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conventional integration method may be used in conjunc- 
tion with a VUV frequency comb [18H20] to compare the 
nuclear clock to an electronic anchor line in the same ion 
or to a 27 A1 + clock, for example, in order to probe for 
variation of fundamental constants at an unprecedented 
level of precision [3T]. With a laser-noise suppressing 
clock interrogation protocol (25], the long excited clock 
state lifetime could be more fully utilized, allowing for 
a pair of 229 Th 3+ clocks to be compared with < 10~ 19 
statistical uncertainty in only a few days of integration. 

In conventional optical atomic clocks, the most signif- 
icant external-field clock shifts arise due to differences 
in the clock state electronic wavefunctions. In the nu- 
clear system proposed here, all four clock state electronic 
wavefunctions are the same. External-field shift analy- 
sis of this compound system (electron cloud + nucleus) 
must then be applied to changes in nuclear properties and 
the resulting slight changes they create in the electronic 
wavefunction. Consequently, the difference for external- 
field shift property X (magnetic-dipole moment, electric- 
quadrupole moment, polarizability, etc..) comes in three 
flavors: 

AX=j x AX nuc + AX tso + AX hfs , (1) 

(i) AX nuc is the direct nuclear contribution including 
core-electron shielding of externally-applied fields, 7x 
being the shielding factor of Th 4+ , (ii) AX tso is the 
isomer nuclear-charge distribution (monopole-mediated) 
shift, and (iii) AX h ^ s is the valence-electron contribution 
to the hyperfine-mediated shift. 

To analyze contributions AX lso and AX h f s , the elec- 
tronic response of 229 Th 3+ to changing nuclear proper- 
ties is computed using techniques of relativistic many- 
body perturbation theory (MBPT). The calculations 
start from the Dirac-Hartree-Fock approximation and in- 
corporate the major correlation effect - core polarization 
by the valence electron. Calculations of matrix elements 
additionally include an all-order chain of random-phase- 
approximation diagrams. It is found that this approach 
recovers experimental hyperfine structure constants with 
an accuracy of 10%, which, considering uncertainties in 
nuclear parameters, is sufficient for the goals of this Let- 
ter. Finally, various shielding factors for the Th 4+ core 
are taken from Ref. [2"3"] . 

To determine the monopole-mediated contribution, 
AX lso , the nuclear charge radius R rm s is varied for 
a series of MBPT calculations of property X. This 
yields the rate of change dX/dR rms . As a result, 
AX lso = dX/dR rms AR rms , where AR rms is the dif- 
ference in the rms radii of the nuclear ground- and 
isomer-state charge distributions. According to nuclear- 
structure calculations [23], \AR rms \ < 0.0038 fm for the 
nuclear transition in 229 Th. 

The hyperfine-mediated shift is determined via first- 
order perturbation theory in the hyperfine interaction of 



electrons with the nuclear magnetic-dipole and electric- 
quadrupole moments. Generally, for the stretched states 

Ax hfs = ( Mm - Vg \ xhfs + ( Qm-Q g \ xhfs (2) 

V J V \ e \ b / 

where /i g , m and <2 g , m are nuclear magnetic-dipole and 
electric-quadrupole moments of the ground (g) and iso- 
mer (to) states, and form-factors X^q depend only on 
the electronic wavefunctions. Electronic form-factors are 
evaluated below using MBPT methods. The nuclear Ml 
moments used are [i g — 0.45 hn [2E) and /i TO ~ -0.08 /Ltjv 
[13] . The difference in intrinsic quadrupole moments of 
the two nuclear charge distributions, AQ20J w &s investi- 
gated in Ref. [HJ. The condition |AQ 2 o| < 0.28 |e|fm 2 is 
taken from these results and, with Q g — 3.11(16) |e|b |17j . 
the isomer quadrupole moment is Q m = 1.74(17) |e|b. 

Zeeman shifts arise due to interaction of exter- 
nal magnetic fields with the atomic magnetic dipole 
moment /i. In a weak magnetic field B, the clock 
transition shifts are AE^=±Afj,B, where m accor- 
dance with Eqn.Q, Afi = 7 A1 A// mc +A/i iso +A^- fs . 
The individual contributions are found to be 
7A1 A/i nuc «-0.53/Mr, |A// so | <8x10~Vjv, and A/i'^w 
-2xlO _3 /iAf, which lead to differential Zeeman shifts of 
AI?^/Bra^4kHz/mT for the two clock transitions. 

Due to the system's symmetry, only a magnetic field 
change which is correlated with alternation between clock 
transitions will produce a clock shift. All other fre- 
quency components, e.g. ac line and /zT-level fields at 
the trap drive frequency, will only modify the clock spec- 
trum with small resolved modulation sidebands or slight 
line broadening. The correlated field change should be 
held to < 10 pT at the ions to maintain a fractional shift 
of lxlO -20 . Magnetic shielding with 10 -3 field suppres- 
sion placed around the vacuum system could be combined 
with microwave spectroscopy, where the magnetic-field 
sensitivity enhancement is /ib/hn ~ 1800, to ensure this 
criterion is met. 

The dominant contribution to second order Zee- 
man shifts is related to the magnetic-dipole polariz- 
ability of the nucleus itself. This quantity is en- 
hanced due to the unusually small energy gap be- 
tween the ground and isomeric states of the nu- 
cleus. As a result, the Ml polarizabilities of the two 

states have the same absolute value but opposite signs. 

(2) 

The relevant clock shift may be estimated as AE Z ' 
= (2\(g\A4 n) \m)\ 2 B 2 )/AE clk ^(2^ N B 2 )/AE clk . A 
1 /iT field fluctuation at the clock ion, e.g. from rf cur- 
rents in the trap electrodes, would then lead to a frac- 
tional clock shift of only ~4xl0~ 29 . 

The quadrupole shift arises due to interaction of the 
atomic quadrupole moment Q with external electric 
field gradients. The value of 7 Q AQ lmc w 9x 10" 6 |e| a% 
is estimated using AQ nuc = -1.37(5) |e|b discussed above 
and the quadrupole shielding factor "fQ — -177.5 for 
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Th 4+ j^Hj. The monopole-mediated shift is calcu- 
lated to be |Ag iso |<8xl(r 7 |e|ag. Electronic form- 
factors entering the hyperfine-mediated shift, Eqn.([2]), 
are found to be Q h J s =-2. Ox 10~ 6 |e| a 2 , and Q h Q fa '=- 
5.1xl0~ 7 \e\al, leading to AQ h ^ s w 1.8xl0 -6 |e| a 2 ,. 
Combining all three sources of the quadrupole clock 
shift, AQ w 1 x 10 -5 |e| Oq, a result dominated by the anti- 
screened direct nuclear contribution. 

The dominant electric field gradient at the clock ion 
is due to the dc electric field which produces axial ion 
confinement within a linear rf trap. Using an axial secular 
frequency of 500 kHz, the resultant 8x 10 6 V/m 2 gradient 
would produce a fractional frequency shift of « 3x 10~ 20 , 
a value sufficiently small to avoid clock averaging over 
three orthogonal quantization axes |26| . 

The Stark shift arises via atomic polarization induced 
by external electric fields. With an exceedingly small 
differential nuclear polarizability (Aa lmc < 10~ 14 Oq) and 
efficient screening of dc electric fields at the nucleus by 
the electron cloud, the direct nuclear contribution to 
static Stark shifts is found to be negligible. Using MBPT, 
the static scalar (s) and tensor (t) electronic polariz- 
abilities are evaluated for the 5F 5 / 2 state, resulting in 
a s s» 13 aj] and at ~ -4.2 a§. From these results with var- 
ied nuclear size, |Aa lso | < 3xl0~ 5 a 3 ,. 

Calculation of the hyperfine-mediated shift, Eqn.Q, 
is done using the formalism and codes in Ref. [37]. The 
resulting form-factors are, in units of a^: [oC ]s—- 
4.9xl0- 6 , [a^ s ]t = 1.4xl0- 5 , [eiQ fs } s = 7.1xl0~ 6 , 
and [qq /s ]( = -9.3xl0~ 6 . With these values, the 
HFI-mediated differential polarizability is Aa h f s w - 
lxlO~ 5 aQ and, finally, |Aa| <4xl0 _5 ao, dominated by 
the electric monopole-mediated contribution. 

The fractional blackbody radiation (BBR) clock 
shift is parameterized as ^VBBRl v cik — /?x(T/T ) , 
where T = 300 K. In atomic units, P=- 
[(a 3 ^ 2 T 4 )/(15^ fe )]x[Aa] 5 ., where [Aa} s is the 
difference in the static scalar polarizabilities of the two 
clock states. From the static polarizability presented 
above, \/3 iso \ < lxl0~ 22 and (3 hfs « 3x 10~ 23 , with the 
direct nuclear contribution being negligibly small. The 
BBR shift is dominated by the electric monopole- 
mediated contribution and translates into a fractional 
clock shift of <1.3xl0 -22 at room temperature, sev- 
eral orders of magnitude below BBR shifts in current 
generation atomic clocks. 

Interrogation with the clock laser will itself induce a 
differential light shift on the clock states. Because the 
nuclear differential polarizability is exceedingly small, a 
direct nuclear contribution is negligible. For a circularly- 
polarized interrogation field, Aa lso (w c jfc) ~ -5.3x 10~ 6 a 3 ) 
within the expected clock transition energy range of 6.1- 
9.1 eV. Similarly, the hyperfine-mediated contribution 
remains nearly constant at Aa h ^ s (u! c ik) ~ -1.6x 10~ 7 a 3 ,, 
leading to Aa(u> c i k ) ~ -5.3x 10~ 6 a 3 ,. 



To evaluate the interrogation laser light shift, a lower 
estimate of 0.2 /in is assumed for the clock transition ma- 
trix element. A cw laser intensity of 10mW/cm 2 would 
then induce a clock Rabi frequency of ~ 1 Hz and result 
in a fractional light shift of ~ lxl0~ 24 . 

A room-temperature vacuum system would house the 
clock ion trap and could be maintained at a pressure 
of lxl0~ lo Pa. With this low background pressure and 
extremely small differential polarizability of the clock, 
background-gas collisional shifts are expected to be neg- 
ligible, though our estimate allows for an uncertainty of 
lxlO" 20 . 

In addition to a narrow clock transition and sizable 
external-field shift suppression, the 229 Th 3+ clock con- 
tains properties favorable for significant reduction of 
second-order motional (time dilation) shifts, a dom- 
inant systematic in ion clocks. The 150 kHz wide 
5-^5/2 6D3/2 electric dipole transition at 1088 nm al- 
lows for direct ground-state cooling to minimize motion 
while thorium's large mass further reduces the time dila- 
tion shift. To maintain motional ground state occupation 
throughout clock interrogation/evolution, a 232 Th 3+ ion 
may be used for sympathetic cooling. The cooling field 
would overlap both ions, cooling the 232 isotope, while 
an additional 1088 nm compensation field would serve to 
mitigate differential light shifts of the clock states. 

A pair of Th 3+ ions would be stored in a linear rf 
trap where strong rf confinement places both ions on 
the trap axis (z-axis). Within a trap of radius 700 /im 
and drive frequency of 25 MHz, the four transverse mo- 
tional eigenfrequencies may be set to 2.0 -2.5 MHz with 
the two axial frequencies at 0.5 -0.9 MHz. In this case, 
sideband cooling well into the ground state could be 
achieved for all modes directly with a 1088 nm field 
[2"51 12"9"] . The total kinetic energy of the clock ion is then 

EKtt^hii/xo+Vxn+Vyo+Vy^ + ^hfrzo+VzT-), where v ig 
is the secular frequency along dimension i with rela- 
tive ion motional phase 0. Time dilation results in a 
fractional clock frequency shift of SvTD/^cik = -Ek /mc 2 ! , 
where m is the ion mass. For Ek associated with 
ground-state secular motion and minimal micromotion, 
5vTD/vdk ~ -5 xl0~ 20 . With motional state measure- 
ments of (riig) <C \ for all modes, the uncertainty of this 
shift could be held to < lxl0~ 20 . 

A second contribution to time dilation comes from ex- 
cess micromotion (EMM) [30] arising from stray elec- 
tric fields which push the clock ion away from the rf 
potential minimum. In the pseudopotential approxi- 
mation, EMM kinetic energy due to a stray electric 
field Ei along radial dimension i may be estimated as 
Eemah ~ (2m)" 1 (q Si /27r^o) 2 , where q is the ion charge. 
If stray electric fields along the x- and y-dimcnsions are 
confined to < 1 V/m with proper trap design and con- 
tinuous motional monitoring, the additional contribu- 
tion to time dilation from excess micromotion would be 
< l.OxlO" 19 . The uncertainty of this systematic shift is 
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taken as the shift upper estimate itself. 

The electric fields used for ion confinement will induce 
a small differential Stark shift on the clock states [3D] ■ For 
the ground state cooled system with a 1 V/m transverse 
electric field, a shift is estimated from the results above 
to be <3xlCT 26 . 

A single a transition at 1088 nm is dipole allowed 
from each of the four clock states. A light shift of the 
i th clock state due to coupling of this transition to the 
232 Th 3+ coo ii ng . nc id i s A^ioss ~ H%/(4Ai), where tt R is 
the resonant Rabi frequency induced between stretched 
states within the cooling ion and A< is the detuning of 
the cooling field from the relevant 1088 nm transition in 
the clock ion. This detuning is dominated by the 10 GHz 
isotope shift [T7j. Though the dipole matrix elements of 
all four 1088 nm transitions in the clock ion are equal, the 
leading differential light shift arises via different values of 
Aj, due primarily to the < 1.4 GHz isomer field shift [31] 
and secondarily to the ~ 100 MHz differential hypcrfine 
shift. The differential shift of the clock transition is then 
< 15% of the nuclear ground clock state light shift. Using 
f2fj = 150 kHz for ground state cooling equates to a clock 
shift of < 85 mHz. 

An additional 1088 nm field can be used for Stark 
compensation, with the leading limitation in suppression 
due to uncorrelated intensity fluctuations between the 
cooling and compensation field. With 10 -3 fractional 
intensity instability between the two fields, differential 
Stark shifts can be suppressed by three orders of magni- 
tude, leading to a clock shift and uncertainty of 5x 10~ 20 . 
Off-resonant scattering from the clock states due to the 
1088 nm cooling and compensation fields is estimated to 
produce < 10 mHz of clock transition broadening. 

To remove linear Doppler shifts which are correlated 
with clock interrogation, alternating clock probes from 
counter-propagating directions may be averaged [TUJ . By 
minimizing changes in optical intensities, polarizations, 
and detunings when alternating between probe direc- 
tions, this differential linear Doppler shift can be held 
to 1 xl0~ 17 [TT]. With an asymmetry of the clock line- 
shapes and respective servo gains held to 0.1%, the frac- 
tional linear Doppler uncertainty can then be confined to 

1X10- 20 . 

Gravitational shifts will play a significant role in clocks 
operating with inaccuracy near lxl0~ 19 . The frac- 
tional frequency difference between two clocks at differing 
heights on Earth's surface is Avq/v^ = g Ah/c 2 , where 
g is the mean gravitational acceleration at the clocks 
and Ah is the clock height difference. This results in 
a lxlO -19 fractional frequency shift for a 1mm discrep- 
ancy in clock height, a likely level of error for two sepa- 
rated clocks in the same lab. 

A large set of clock shifts relevant to the proposed sys- 
tem has been analyzed, with an error budget including 
these shifts and uncertainties found in Table HI With all 
classes of shift uncertainty added in quadrature, the total 



TABLE I. Estimated systematic error budget for a 229 Th 3+ 
clock using realized single-ion clock technologies. Shifts 
and uncertainties are in fractional frequency units (Av/v c ik) 
where v c ik — 1.8 PHz. See text for discussion. 

Effect [Shift! (1CT 20 ) Uncertainty (ltr 20 ) 



Excess micromotion 


1U 




Gravitational 





10 


Cooling laser Stark 





5 


Electric quadrupole 


3 


3 


Secular motion 


5 


1 


Linear Doppler 





1 


Linear Zeeman 





1 


Background collisions 





1 


Blackbody radiation 


0.013 


0.013 


Clock laser Stark 





<0.01 


Trapping field Stark 





<0.01 


Quadratic Zeeman 








Total 


18 
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estimated clock inaccuracy is 1.5x10 1 . 

Utilization of the 229 Th nuclear transition for high- 
accuracy frequency metrology would require ultra-stable 
interrogation light at 163(10) nm. With a //W- level cw 
power requirement, such light might be generated by 
multi-stage sum-frequency generation of visible or NIR 
light within nonlinear crystals, e.g. LBO, BBO, and 
KBBF, with the fundamental frequency locked to an 
ultra-stable optical resonator. To determine the state of 
the clock after interrogation, a double-resonance method 
may be used where NIR or visible fluorescence from an 
electric-dipole transition in the clock ion is correlated 
with a particular nuclear state [T^] . Due to the long iso- 
mer state lifetime, the clock interrogation protocol could 
begin from both the nuclear ground (excitation) and iso- 
mer (de-excitation) manifolds to expedite the clock ini- 
tialization process. 

In conclusion, this work has shown that the usual 
sources of clock shifts are suppressed in the proposed sys- 
tem, so that clock operation with inaccuracy approaching 
lxlO -19 appears viable. With a more aggressive pur- 
suit of systematic shift compensation and suppression, 
this value may very well be lowered to the 10 -20 scale, 
however, other systematic clock shift mechanisms may 
become operational at this level. 

This work was supported by the Office of Naval Re- 
search, the National Science Foundation, and the Gordon 
Godfrey fellowship. 
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